Squeezed coherent state undergoing a continuous 
nondemolition observation 



J2' 



Anita Dabrowska, Przemyslaw Staszewski 

Department of Theoretical Foundations of Biomedical Sciences and Medical Informatics 
Collegium Medicum, Nicolaus Copernicus University 
ul. Jagiellohska 15, 85-067 Bydgoszcz, Poland 



Abstract 

The time evolution of a squeezed coherent state conditioned by the results of a single 
and double heterodyne measurement is discussed. The mean values of quadratures as 
well as the dynamics of quadrature uncertainties have been obtained within the frame- 
work of the theory of continuous measurements based on filtration equations. It has 
been found that while the mean values depend on the measured noise, the uncertain- 
ties in the optical quadratures are deterministic. Explicit solutions for the latter have 
been provided. Finally, a time development of the squeeze parameter for the posterior 
squeezed coherent state has been found. 
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1. Introduction 

The problem of continuous measurements in quantum systems is one of the most 
challenging fundamental issues of modern theoretical physics |lLl2l]- The quantum fil- 
tering theory developed by Belavkin JsKH makes it possible to describe the dynamics 
of a quantum system continuously observed in time. The role of a measuring apparatus 
is played here by a Bose field which in quantum optics can be treated as an approx- 
imation to the electromagnetic field. The interaction between the quantum system in 
question (system S) and the reservoir (the Bose field) is taken in the Markovian approx- 
imation, i.e. the correlation time of the reservoir is much shorter than the time scale of 
the dynamics of S. The continuous trajectory of results of the observation of the output 
Bose field (Bose field after interaction with S) determines the conditioned evolution of 
S (the posterior state). The main motivation for determining the conditional state of a 
quantum system is developing methods of the quantum feedback control. The quantum 
trajectory has been essential to the design of a quantum control algorithm @|. 
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In this article, we present the posterior evolution of an optical cavity mode coupled 
to the outside radiation mode by one or two partially transmitting mirrors. The time- 
development of a state of an indirectly observed system is conditioned by a trajectory 
of the results of a single and double heterodyne measurement. We assume the coupling 
system's operator proportional to the annihilation operator. In the paper we discuss the 
analytical solutions to the filtering equation for the system being initially in a squeezed 
coherent state. In order to prove that such a state is preserved under the considered 
continuous diffusion observation we use the linear version of the quantum filtering 
equation derived in yj. 

The exact solutions to the filtering equation for the initial Gaussian states and a dif- 
fusion observation were given, for instance, in I7I- H0I1 . The solution to the quantum 
filtering equation for a harmonic oscillator in an arbitrary initial state state undergoing 
the heterodyne observation was discussed by Carmichael in 111 111 . 

The main part of our work is organized as follows. In Section 2 the mathematical 
model is described. Section 3 contains the analysis of the double heterodyne detection 
problem. Section 4 is devoted to the single heterodyne detection. Some final remarks 
are given in Section 5. 



2. Model 

We consider a single cavity mode of the electromagnetic field (a system S) interact- 
ing with two independent components B„(t) , (n — 1, 2) of the Bose field being initially 
in the vacuum state. The unitary operator, U(t), describing the evolution of the whole 
system (the system S plus Bose field) satisfies the Ito quantum stochastic differential 
equation (QSDE): 



dU(t) = 



h ^4 



n=\ 



At 



+ Z yj^{adBl(t)-aUB n (t) 



U(0) = / 



U(t), 



(1) 



where H - ha> (a^a + is the hamiltonian of S, a stands for the annihilation operator 
of S and fi„ > (n = 1, 2) are coupling constants. Eq. (Q} is written in the interaction 
picture with respect to the free dynamics of the Bose field. The discussion of the 



physical assumptions leading to this evolution one can find, for instance, in [ 121 11311 . 
In brief, the interaction hamiltonian is taken linear in the field operators, the rotating- 
wave approximation (RWA) is made and a flat and broad spectrum of the reservoir is 
assumed. 

The Bose field provides a possibility of a continuous indirect observation of S. We 
suppose that the information about the systems is gained by using a double heterodyne 
detection scheme, thus we consider a simultaneous measurement of the two output 



2 



processes 



i 

o 

+ ^[^IJ2Rt[t-' iM ' ) a)At' , n=l,2, (2) 

where cf>(t) — (pQ + §t Il4l[l5ll . These output processes satisfy the Belavkin's nondemo- 
lition condition, namely 

[Q™\s),U~Ht)ZU{t)] = Vs<f, (3) 

where Z is any operator of the system S. Furthermore, the output processes (|2), due to 
their Hermicity and self-commutativity, 

[Q° ut (t),QT(t')] = V?,f'>0, (4) 

can be treated as the classical Wiener processes. 

The time development of the posterior unnormalized wave function t//(t) of S, corre- 
sponding to the trajectory of the observed processes (O up to t, is given by the Belavkin 
linear filtering equation of the form 

2 



di//(t) = -fr# + Z ^a f a\i//(t)dt 

\ n=l ' 



+ Z V T a e ~ lM ° m dQn(t) ' ^ (0) = ^' (5) 

n=\ " 2 

where <2„(f) (n = 1,2) are Wiener processes for which d<2„(f)d<2„,(f) = 6„„,dt. Let us 
recall that the posterior mean value of any operator Z of the system S reads in terms of 
unnormalized wave function if/it) satisfying the linear filtering equation 

<Z>, = — — — . (6) 

WMt)) 

For more details on the role of the linear filtering equation for the diffusion observation 
in quantum mechanics, its derivation, and some exemplary solutions see for instance 



in quan 



3. Double heterodyne detection 

Let us put: fi„ = fi (n — 1,2), 02(0 - 0i(O = jr /2, and <pi(t) = <p(f). With these 
assumptions, the linear filtering equation for the double heterodyne measurement takes 
the form 

d4(t) = ~{^H+ ja t aj/f(f)dt + y/p a e~ m ^(t) dQ(t) , 

m = f, (7) 
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where d<3(f) = -^= (dQ\(t) — id(32(0) is the complex diffusion process such that 
2 and|d<3 + (f)l 2 = dr. 



(dfl(0) 



Let us discuss the time development of the posterior wave function, assuming that 
the initial state of S is a squeezed coherent state [19], 



where 
and 



<A(0) = S({o)D(a )\0) = S(ft)|a > = lft,<*o>, (8) 
D(ao) = exp (aoa ' - aofl) , (9) 

S(ft) = exp^i 0fl 2 - j , ft = e lflo eo e C. (10) 

The amount of squeeze is described by the modulus po of the squeeze parameter go, 
whereas the phase Gq specifies the angle of the squeeze axis in the phase space. Ob- 
viously, for £o = the state (O becomes the coherent state \ao), and for ao = and 
ft + it takes the form of the squeezed vacuum state, S'(ft)IO). Let us recall that the 
expectation values of the quadratures vanish while the expectation value of the photon 
number operator is nonzero in the state 5(ft)|0). 

We shall prove that the solution to Eg. (O corresponding to the initial state © can 
be written as 

m = i(t)s(mw)). (id 

The proof is rather standard but still quite cumbersome, therefore we shall present its 
outline augmented with some computational details and partial results. 

Substituting the postulated solution (fTTT > into Eg. (O, writing both sides of the equa- 
tion in terms of linearly independent vectors J la), } and comparing the coeffi- 
cients of the corresponding vectors yield the consistent system of differential equations 
for the functions: g(t), 6(t), a(t), l(t). The solution to this system corresponding to the 
initial state ([8]l uniquely defines the solution to Eg. (0 in the form ( fTTT ). Let us pay 
attention to several steps of these calculations. 

To calculate the increment di/r(f) one has to calculate the increment dS (£(?)) = 
5 (£(/ + df)) - S 1 To this end it is convenient to make use of the normally ordered 
form i20ll of the squeeze operator, S (£), 

S(£) = (cosh q)~ 1/2 exp -r(a f ) 2 /2 

x exp [-ln^osh^fl'fljexpfra 2 ^] , (12) 

where 

r = e ifl tanh£). (13) 
Left-multiplying Eq. (0 by S ' (£(0) and using the unitary transformation 

sU&^'fs = (a t ) 2 cosh 2 g 
- (2a"'a + l)e~ ifl sinhecoshe + aV 2ifl sinh 2 p, (14) 
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one gets 



S\€(f))dS(£(t)) = i tanhg(t)dg(t) [2T(t)a 2 - 2a 1 a - l] 



H — - — a - 



2 

dT(t) 



cosh 2 p(f) 



2 2 

- (2a + a + 1) e~ ifl(,) sinh £.(0 cosh g(t) +a V 2ifl(0 sinh 2 £>(?)] . 
By the virtue of applying the formulae 



= a cosh q - a ' e 1 sinh^i, 

= fl'acosh 2 £) + (a' a + l)sinh 2 £ 



a 2 e- ie ■ 



(4 



sinh £> cosh p 



and 



a' \a) 

(«"'")V> 



d\a) 1 
^ + 2 a|a> ' 
d 2 |a) | _g|g) | 1_ 2 
5a 2 3a 4 



(15) 

(16) 
(17) 

(18) 
(19) 



one is able to rewrite the equation in terms of the vectors ^r}. As these 

vectors are linearly independent, the comparison of the expansion coefficients yields 
the set of differential equations for the functions 0, g, a, and I: 



d6(f) 
da(t) 

dl(t) 



= -2a)dt , dg(t) = -fi sinh g(t) cosh g(t)dt , 



- \ia) + ~ j - ju sinh 2 g(t) 



a(t)dt 



- V^e lfl(,) sinhe(f)e" lW, d(3(f) , 
-ydf + Id \a(t)\ 2 + £ sinh 2 g(t) \a(t)\ 2 dt 

+ ^7ff(f)cosh £) (f)e" iw d<3(f) - | sinh 2 £>(0df 

+|a 2 (f)e _iS( ' ) sinh£>(0 coshe(f)df , 

with the initial condition 1(0) = 1, a(0) = or , #(0) = #o, £>(0) = go- 
The solution to the system (EOt reads 



(20) 



0(f) 
o(t) 



= e 



0o - 2u)t , g(t) = artanh (e A " tanhpo) 
(teH-S)i coshg(0/_ 



- , 

coshpo \ 

f 

V^e ie » sinhpo J e-^)V W) d<a(0 



(21) 



/(r) 



cosh p(f ) 



V cosh^o 



■exp 



^ + U\a(t)\ 2 -\a \ 2 )+ X (t) 
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where 



x (t) = yfJlJ{a(t')coshQ(t')e- W) dQ(t') 
o 

+ -y£<e- ifl( 'V(f') smhg(t') coskg(f)dt') . 

Therefore the posterior mean values of optical quadratures X — (a + a T )/2 and Y = 
(a - a T )/2i for the posterior wave function of the form (fTTT i, given by 

(X), = Re (a (r) cosh g(t)-a (?) e ie(0 sinh g(t)) , (22) 

(Y), = Im (a(t ) cosh g(t) + a(t)e-' ,m sinh g(t )) , (23) 

depend on the measured noise, whereas the uncertainties in the optical quadratures 
AX(t) and AT(f) are deterministic. One has: 

AX(t) = ~ [l + C(0 (e^' tanhpo - cos (0 O - 2wf))]' /2 , (24) 
AF(0 = i [l + C(f) (e^' tanh eo + cos (6» - 2wf))]' /2 , (25) 

where 

_ 2e"^ tanhgo 
1 - e- 2 ^' tanh 2 go ' 

The squeezing of X occurs for 

cos (6»o - 2a>t) > eT>" tanh go , (26) 
while the squeezing of Y occurs for 

cos (0q - 2a)t) < -e'"' tanh g . (27) 

The time dependence of the uncertainties AX and AT has been illustrated by the 
parametric plots presented in Fig. 1. They show the dynamics of AX and AT as func- 
tions of the dimensionless time r = tot (0 < r < 100) for /i = O.Olw #o = 0, and 
three values of po. The pictures for different values of po differ mostly for short times, 
otherwise they are qualitatively very similar. From the above formulae as well as from 
the figure it is clear that the dynamics of quadrature uncertainties is such that the sys- 
tem switches back and forth from being squeezed in one of the quadratures but not 
in the other while passing also through the region where there is no squeezing at all. 
Asymptotically, the system approaches the vacuum state (with AX = AY = i). An in- 
teresting feature of the envelopes of the displayed curves is that the region where they 
are concave shrinks as time grows, and the envelopes finally become convex. 

A similar result concerning the uncertainties in the posterior optical quadratures as 
well as the squeezing coefficients was obtained in [16] where subharmonic generation 
from the vacuum was studied: the uncertainties and 77(f) were found to be deterministic, 
while a(t) was generated from the output noise. 
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Figure 1 : Time dependence of the uncertainties AX and AY as given by Eqs. d24!25t . The dependence of AX 
and A7 on the dimensionless time t = ait is displayed {or /i = O.Olu, do = 0, & = 0.05 and for three values 
of p : 0.5 (a), 2.0 (b), and 8.0 (c). 
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It may be appropriate at this point to recall that the posterior wave function satisfy- 
ing the linear version of the filtering equation is normalized to the probability density 
of the output (observed) diffusion process with respect to the standard Wiener measure 
of the input Wiener diffusion process. In the considered case of the double heterodyne 
observation Eq. (fTTb implies ||i^(f)ll 2 = | /( ? ) | 2 , therefore \l(t)\ 2 is the probability density 
of the complex Wiener process Q(f) with respect to the standard Wiener measure of the 
complex input process Q(t). Note that the function 1(f) is essential in the considered 
continuous observation. For example, as one can see from d2Tb that l(t) depends on all 
the parameters of an initial squeezed coherent state. If the values of some of them or 
even all the values of these parameters are unknown, they can be determined form l{t). 

We have proved that the filtering equation makes it possible to study the time- 
development of a squeezed coherent state. This would not be possible with the help 
of the master equation which can be obtained from the filtering equation for the mixed 
posterior state obtained from the Eq. (0 by taking the stochastic average over all pos- 
sible trajectories of the observed process. Though the master equation (of the same 
form for both cases of the heterodyne observation considered in the paper) preserves 
an initially coherent state [ 17], it does not preserve the squeezed coherent one [18]. 

4. Single balance heterodyne detection 

The filtering equation for a single balance heterodyne detection 

dfi(t) = -[~H + ^a'flj^Odr + y/ji a e" w ^(r) d<3(0 , 

MP) = iff, (28) 

can be obtained from Eq. (0 by putting /j\ = 2/j., /j 2 = 0, <p\(f) — <p(t), Q\(t) = Q(t). 

We shall prove that the squeezed coherent state is preserved under the diffusion 
observation. For this purpose we employ the property 

S(g)aS\&\S,a) = a\t,a). (29) 

Making use of the Baker-Hausdorff formula one can get the unitary transform of the 
operator a 112011 

S(£)aS\g) = aTi+cJTi, (30) 

where Fi = coshg, T 2 = e lfl sinhp. 

Let us notice that if the system remains in the squeezed coherent state ( fTTb at any 
time instant t > 0, then the following relations have to be satisfied 

sm))as\mm = [ fl r 1 (f) + « t r 2 (f)]^(f), od 

S (£(t + dt))aS f (£(t + dr)>A(f + At) 

= [a Fj (f + df) + a f T 2 (t + df)] ty{t + df) . (32) 
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Eqs. (|3TT > and d32l can be reduced to the single condition 

[a (Ti(/)+dr 1 (*))+a t (r 2 (0+dT 2 (0)-t!r(0~dar(0] d?(0 
+ (adr 1 (f) + a + dr 2 (0-da(f))f(f) = 0. 



(33) 



Then by insertion of di/^f) from Eq. d28b into Eq. d33l we obtain the set of the differen- 
tial equations 



o(t)ri(0 



-r, (f) d/+dTi(0 +ju e _2W r 2 (0d/ 



-dar(0-a(0r 2 (0 



2 (o(iw+|WdT 2 (0 



(34) 



V^r 2 (f)e^ w d(3(f) = 



T 2 (f) 

-r,(o 



-ri(f) (ioj + -\dt + dTi(f) +//e _2W r 2 (0df 

(t)(k» + |)d* + dr 2 (0 



(35) 



with the initial condition: F^O) = coshpo* r 2 (0) = e' e °sinh£)o> <*(()) = ao, and this 
completes the proof. 

The equation d35b imply that the function F(f) = e ie ® tanhp(f) satisfies the Riccati 
differential equation of the form 



lr(o = -2(i w + ^)r( f) + //e- 2 ^>r 2 ( f ), 



T(0) = e ie °tanh£>o. 



(36) 



Hence the posterior uncertainties of quadratures for the posterior squeezed coherent 
state, given by the formulae 



AX(t) = (4Re*(f)r 1/2 , 
AY(t) = \K(t)\ (4Re*(f)r 1/2 



where 



1 + iXO 



(37) 
(38) 

(39) 



1-1X0 

do not depend on the measured noise, as before. The general solution to Eq. d36l l can 
be written as 

r(0)e- (2iw+/j)f 



i-^r(0)/ei 2iw+ ^''^ ( '' ) df 

o 

In particular, for the phase 0(f) = n/2 + §t, we obtain 

(2iw + 2i§ + fx)T(0) 



(40) 



no = 



e( 2i ^'"[2iw + 2i# +//(! + r(0))] - y. T(0)e 



-2i#( 



(41) 



w 111 1 



1 1.5 2 2.5 3 3.5 




AX(t) 



Figure 2: Same as Fig. 1 except for Eqs. (36-37) with <p(f) = n/2 + -&t. 



The time dependence of the uncertainties AX and AY computed above has been 
illustrated in Fig. 2 showing the dynamics of AX and AY (with the same values of all 
parameters as in Section 3). It is very hard to find any qualitative difference between 
the corresponding parts of Fig. 2 and Fig. 1. Such a difference exists only for small 
times. This is because the crucial variable F(f) becomes asymptotically exponential for 
large times, that is, it behaves in the same way as in the case considered in Section 3. 

To derive the differential equation for the coefficient l(t) one has to insert the state 
l(t)\^(t), a(t)) to the filtering equation ( f28l ). The calculations are more complicated in 
comparison with these of Section 3. Finally one gets: 



/(f) = exp - ^ + l -{\a{t)\ 2 - |a () | 2 ) 



r'( moi 2 nt)a 2 (t) " 
Jo U(i-|r(0i 2 ) i-moi 2 , 



df 
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with 



H f ' (r(0) 2 e 2i ^ w / 1 



2 Jo 1 

f 

Jo 



(42) 



o(0 



>/l - W(t)\' 



dQ(t) 



a(t) = 



i - HW 



fl'O 



mo)| 2 



x exp -lwf 1 + n 



-20(f) 



r(?')df 



Jo 

^I exp (-Hf) 



(t-s) 



(43) 



5. Final remarks 

We have shown that in contrast to the master equation, the filtering equation de- 
scribing the reduction of the quantum state following the registered trajectory for a 
single and double heterodyne measurement does not destroy the squeezed coherent 
state (fTTt . We have checked that for the initial state |£o, ao) with £o + 0, the amount 
of squeezing decreases in time and registering the trajectory increases our knowledge 
about S. It is also worth to emphasize that for the system prepared initially in a co- 
herent state, the prior and posterior mean values of the system operators coincide. The 
posterior coherent state has a random phase and decreasing in time independent of the 
noise amplitude. Even this case one can take advantage from using the filtering equa- 
tion: the probability density \l(t)\ 2 depends on ao, so l(t) gives information on the initial 
state of the oscillator. The analytical solutions to the filtering equation for the initial 
coherent state obtained for the diffusion observation in 3 21], and for the count- 
ing process in 12211 are consistent with our results. Finally, let us note that it will be 
very interesting to consider the posterior evolution of a squeezed coherent state when 
a driving force is applied to the oscillator. 
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